Large magnetothermal conductivity in GdBaCo 2 05 +x single crystals 
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To study the effects of paramagnetic spins on phonons, both the in-plane and the c-axis heat 
transport of GdBaC^Os+ic (GBCO) single crystals are measured at low temperature down to 0.36 
K and in magnetic field up to 16 T. It is found that the phonon heat transport is very strongly 
affected by the magnetic field and nearly 5 times increase of the thermal conductivity in several 
Tesla field is observed at 0.36 K. It appears that phonons are resonantly scattered by paramagnetic 
spins in zero field and the application of magnetic field removes such strong scattering, but the 
detailed mechanism is to be elucidated. 

PACS numbers: 66.70.+f, 72.20.-i 
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I. INTRODUCTION 

In magnetic materials, coupling of spins with other el- 
ementary excitations is a matter of fundamental impor- 
tance. In particular, the spin-phonon interaction causing 
unusual physical properties in correlated oxides has been 
an issue of significant interest both by its own right and 
because it provides a useful means to gain insight into 
frustrated magnetism as well as multiferroicsi 1 ' 2 ' 3 ' ' The 
low-temperature thermal conductivity is one of the most 
effective probes for studying the spin-phonon coupling, 6 
especially in insulators where magnetic excitations along 
with phonons are the main heat carriers. Even if the long- 
range magnetic order is not established in a crystal or 
the contribution of magnons to the thermal conductivity 
is limited at low temperatures, the presence of localized 
spins still can play an important role in the heat trans- 
port as was observed long time ago in some ionic com- 
pounds doped with magnetic impurities^^ However, in 
strongly correlated materials, the low-T behavior of ther- 
mal conductivity can be so complicated that it is difficult 
to selectively study the effect of spin-phonon interaction, 
especially when the electron transport, which can also 
depend on magnetic field, is contributing to the thermal 
conductivity. In order to gain insight into the impact of 
spin-phonon coupling and its underlying mechanism, one 
needs a material in which the spin-phonon coupling can 
be treated independently. 

For this purpose, recently synthesized layered cobaltite 
GdBaCo205+a; (GBCO), which has attracted a great 
deal of attention due to its rich phase diagram and fas- 
cinating physical properties (See Ref. and references 
therein), would be an ideal system for studying the spin- 
phonon coupling. This compound is a good electrical 
insulator at low temperatures for a wide range of oxy- 
gen content, and hence there is no need to consider the 
electronic heat conduction. Although GBCO shows a 



long-range magnetic ordering of Co spins, their Ising-like 
spin anisotropy prevents the low-energy magnon excita- 
tions, while the Gd ions remain paramagnetic and be- 
have as almost free spins down to very low temperatures. 
Furthermore, a lattice disorder in GBCO can be easily 
introduced by merely changing its oxygen content. 

In this work, we perform detailed studies of the low- 
T heat transport of GBCO single crystals with various 
oxygen contents. It is found that both the in-plane and 
the c-axis thermal conductivities (K a b and k c ) show ex- 
tremely strong magnetic-field dependence at low temper- 
atures. The k(H) isotherms at 1-2 K show a dip-like fea- 
ture at low fields where the heat conduction is reduced 
by as much as 50%, while at subkelvin temperatures the 
magnetic field induces a step-like enhancement of k by 
nearly 5 times. The magnetic-field dependences of the 
thermal conductivity indicate that phonons are strongly 
scattered by the Gd spins that are susceptible to the Zee- 
man effect in magnetic fields. The large magnetothermal 
effect illuminates the important role of free spins in the 
heat transport of novel oxide materials. 



II. EXPERIMENTS 

High-quality GdBaCo205 +:E single crystals are grown 
using a floating-zone technique^ ' 10 : 11 and carefully an- 
nealed to tune their oxygen contents to x = 0.50, 0.30 
and 0.00. The thermal conductivity along both the ab 
plane and the c axis is measured from 0.36 to 5 K in 
a 3 He refrigerator by using a "one heater, two ther- 
mometers" technique and from 5 to 300 K in a 4 He 
cryostat by using the Chromcl/Constantan thermocou- 
ple. Details of the measurements of the temperature and 
magnetic-field dependence of k were described in pre- 
vious publications ^ 2 : 13 ' 14 All the parallelepiped crystals 
are precisely cut along three crystallographic axes with 
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FIG. 1: (color online) The low-temperature specific heat of 
GdBaCo2C>5.oo in magnetic fields of 0, 5, 10, and 13.5 T ap- 
plied along the ab plane. Solid lines show calculated Schottky 
contributions of paramagnetic Gd 3+ ions assuming S = 7/2 
(see text). The dashed line indicates the Schottky contribu- 
tion of paramagnetic ions with S — 1/2 at H — 13.5 T. The 
inset shows the heat capacity in the temperature range of 2- 
300 K measured at H = and 10 T. The arrow indicates the 
charge ordering transition. 



an error of less than 1°. The heat capacity is measured 
by the relaxation method in the temperature range from 
2 to 300 K using a commercial physical properties mea- 
surement system (PPMS, Quantum Design). 



III. RESULTS AND DISCUSSIONS 

Before presenting the heat transport results of GBCO 
crystals, we show in Fig. 1 the representative heat capac- 
ity data, which prove to be very informative regarding the 
magnetic excitations. In GBCO, regardless of the oxy- 
gen content, the spins of cobalt ions undergo magnetic 
ordering at rather high temperatures-*^ but Gd ions, in 
contrast, remain paramagnetic down to very low temper- 
atures as suggested by magnetization measurements^ In 
order to confirm this, we have measured the heat capacity 
of an x — 0.00 single crystal in different magnetic fields. 
As can be seen in the inset of Fig. 1, the C P (T) curve 
is insensitive to the magnetic field, except for the low-T 
part: the peak around 248 K, which corresponds to the 
charge ordering transition at this temperature^ is not 
affected by applying a high magnetic field either. At low 
temperatures, on the other hand, the specific heat C p /T 
shows a pronounced Schottky anomaly with a maximum 
that shifts to a higher temperature with increasing mag- 
netic field. All the data measured in H = 0, 5, 10, and 
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FIG. 2: (color online) Temperature dependences of the ther- 
mal conductivity along the ab plane and the c axis in zero 
magnetic field for GdBaC^Os+ic single crystals with various 
oxygen content. The data are displayed in both linear and 
logarithmic plots. The dashed lines indicate the T 3 depen- 
dence. 



13.5 T are well fitted by a simple Schottky contribution of 
paramagnetic Gd 3+ ions with a spin value of 7/2, which 
is determined by the area under the peak (the Schottky 
contribution of paramagnetic ions with S — 1/2 at H = 
13.5 T is shown for comparison), and the g factor of 1.64 
±0.04, which is determined by the peak position. We 
should note also that there is a finite splitting of the Gd 
spin levels even in the absence of an external magnetic 
field that can be a result of a weak magnetic interaction 
of Gd ions with each other or with the cobalt sublattices. 

Figure 2 shows the temperature dependences of K a {, 
and n c of GdBaCo205+2; single crystals in zero field. 
First of all, as we already pointed out, the insulating 
ground state and the Ising-like spin anisotropy of GBCO 
leave little room for carriers other than phonons to con- 
tribute to the heat transport, at least at low tempera- 
tures. As can be seen in Fig. 2, the x = 0.00 crystals 
show a large phonon peak at ~ 20 K, indicating the neg- 
ligible imperfection of the crystal lattice. In contrast, 
for x — 0.30 and 0.50, the phonon conductivity is dra- 
matically weakened at low temperatures and the phonon 
peak is wiped out almost completely, which is not sur- 
prising since the oxygen non-stoichiometry brings strong 
lattice disorder. Between the latter two compositions, 
the x = 0.50 crystals have better heat conductivity, pre- 
sumably due to the oxygen ions being ordered into chains 
at this oxygen contenti*^ The worst heat conduction is 
observed in x — 0.30 crystals, which have the largest 
degree of oxygen-induced lattice disorder. 9 
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The temperature dependences of k in Fig. 2 are ac- 
tually more complicated than those in conventional in- 
sulators. At high temperatures, for example, the k(T) 
data for the good phonon conductor x = 0.00 cannot be 
described by the usual phonon Umklapp scattering; 6 -^ 
which suggests the existence of a more complicated scat- 
tering mechanism or some other kind of heat carriers. 
In fact, the specific heat data in Fig. 1 indicate that 
above ^100 K, magnon excitations of the Co spin lat- 
tice become important and may affect the high-T heat 
transport properties. On the opposite side, the lowest-T 
behavior of the heat conduction also seems to differ from 
the boundary scattering limit of phonons expected for 
conventional insulators^ in which case the magnitude of 
the phonon thermal conductivity should depend only on 
the sample size, and not on impurities, defects, or oxygen 
concentration (the dispersion of acoustic phonons and the 
sound velocity do not change much upon changing the 
oxygen content). Although both n a t> and k c below 1 K 
are rather close to ~ T 3 , the heat conductivity of the x — 
0.00 sample is nearly one order of magnitude larger than 
that of x = 0.30 and 0.50 samples. Note that the k c data 
for different oxygen contents were collected on one and 
the same sample that was annealed each time to obtain 
the required x; the samples for K a & differ in size by less 
than 20%. In addition, the low-T parts of the k(T) data 
are not smooth: regardless of the oxygen content, all the 
crystals exhibit clear wiggles in their k(T) curves below 2 
K, which suggests the existence of some kind of resonant 
phonon scattering. 6 

Upon applying the magnetic field at low temperatures, 
we have found remarkably strong changes in the heat con- 
ductivity of all crystals. For instance, at the lowest tem- 
perature 0.36 K, both K a b and k c of the parent compound 
x = 0.50 show a step-like enhancement by several times 
with increasing field and saturates above ~5 T [Figs. 3(a- 
d)]. At higher temperatures, a dip- like feature appears 
at low field, while the subsequent enhancement at higher 
field gradually weakens with increasing temperature. At 
the "high" temperature of 5 K, the magnetic-field de- 
pendence is already rather weak and only a broad and 
shallow dip remains. Here, the most impressive result is 
the magnitude of the field-induced changes in k: the low- 
field suppression can be as large as a factor of 2, and the 
high-field enhancement can be nearly 400%. While the 
magnitude of the high-field enhancement in k depends 
on the direction of the heat current [Figs. 3(a-f)], it is 
surprisingly insensitive to the field direction, which con- 
firms that magnon excitations in the long-range ordered 
spins cannot be relevant to the low-T heat transport. 

A clue given by the dip- like feature in k(H) is that the 
low-T heat transport may be governed by the phonon 
scattering by free spins jSiLS as we previously observed in 
Pri. 3 Lao.7Cu0 4 (PLCO)^ Qualitatively, the k(H) be- 
havior shown in Fig. 3 can be well understood in terms 
of magnetic scattering of phonons: (i) In zero field, the 
phonons are presumably scattered by the spins of Gd ions 
that behave like S = 7/2 free spins judging from the low- 
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FIG. 3: (color online) (a-d) Magnetic-field dependences of 
the low-T thermal conductivity of the GdBaCo2 05+i, single 
crystal with x — 0.50. Directions of the heat current (Jh) 
and the magnetic field (H) are shown in each panel. (e,f) 
Temperature dependences of K a t and k c in zero and 16 T 
magnetic fields. The dashed lines indicate the T 3 and T 
dependences. 



T magnetization and specific heat data^^ To be capable 
of scattering phonons, the Gd spins should not be com- 
pletely free, but the spin states should be slightly splitted, 
which has been confirmed by the heat capacity data, (ii) 
In magnetic fields, the energy splitting of the spin states 
is increased by the Zeeman effect. The phonon scattering 
off these spins is most effective in suppressing the heat 
transport when the energy splitting is equal to ~ 3.8/csT, 
where the phonon conductivity spectrum (defined below) 
peaks; therefore, the spin-phonon scattering generates a 
dip- like feature in n(H) at this energy and the dip posi- 
tion shifts to higher fields with increasing temperature. 6 
(iii) In the high-field limit, the spin energy splitting be- 
comes too large to exchange energy with phonons and 
the spin-phonon scattering is completely quenched, en- 
hancing k above its zero-field value. 

To illustrate the above picture, in Fig. 4 we show 
schematically the relation between the energy splitting of 
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FIG. 4: (color online) (a) Schematic picture of the phonon 
conductivity spectrum n(u>) at several temperatures (the 
curves are normalized by their peak values and shifted verti- 
cally for clarity). The phonons can be strongly scattered by 
magnetic ions if their energy is close to the energy splitting 
of the spin states, which gives rise to a "resonant scattering 
band" centered at A = Ao + gfiBH, where Ao is the zero- 
field splitting and gfJ-BH is the Zeeman splitting, (b) Simu- 
lated k(H) by assuming that the phonons within the resonant 
scattering band do not carry heat. 
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FIG. 5: (color online) Representative data for the magnetic- 
field dependences of the low-T thermal conductivity of 
GdBaCo2 05+:E single crystals with x = 0.00 and 0.30. 



spin states and the peak in the phonon conductivity spec- 
trum from the Debye model. The "phonon conductivity 
spectrum" k(lo) represents the contribution to thermal 
conductivity from the phonons with an energy huj and is 
proportional to w 4 exp(^/fc B r)/T 2 [exp(^w/fc B T) - l] 2 
which peaks at energy hui ~ 3.8/csT. 6 Although Gd ions 
have the spin state S = 7/2 that splits into eight energy 
states in magnetic fields, allowing multiple transitions, 
for simplicity we take just a single band, corresponding 
to the spin transition with the lowest energy, Am = 1. 
The energy of this transition changes with the magnetic 
field as A = Ao + gj-isH, and the phonons in a certain 
energy range around this A are considered to be reso- 
nantly scattered upon the spin transition. Figure 3(b) 
shows simulated n{H) behaviors, where one can see the 
variations of n as the position of the resonant scattering 
band is changed by the magnetic field. Upon calculating 
these k(H), we simply assume that those phonons whose 
energy lies in a resonant scattering band centered around 
A cannot contribute to the heat transport, and the rest 
of the phonons have the same relaxation time (which is 
the case at very low temperature) It is reassuring that 
such a crude picture with only two adjustable parame- 
ters, Ao (~2 K) and the width of the resonant band, can 
capture most of the qualitative experimental features of 
k(H). Indeed, both the step- like enhancement and the 
dip-like feature are essentially reproduced. 

Nonetheless, there is an obvious difficulty in the above 
simple picture; namely, the magnitude of the observed 
changes in k(H) requires the resonance scattering band 
to be extremely wide; indeed, the observed five-fold in- 
crease in k implies that at least 80% of phonon conductiv- 
ity spectrum are affected by the spin scattering. Appar- 



ently, including other resonance bands related to higher 
energy spin transitions of Gd ions with Am > 1 does 
not help, since the Zeeman splitting Amg/isH should 
quench these transitions at several times lower fields than 
we observe in the experiment. A possible explanation for 
the extremely wide scattering band is the local fluctua- 
tion of the zero-field splitting Ao , which would be quite 
natural if Ao is determined by the dipole-dipole interac- 
tions within the disordered Gd subsystem. Depending on 
the local environment, Gd ions feel different local fields 
and thus should exhibit a continuous spectrum of spin- 
transition energies, rather than a single resonance line* 1 ^ 

We have done the same experiments oni = 0.30 and 
0.00 crystals and show the representative n(H) isotherms 
in Fig. 5. Obviously, the main features, including both 
the high-field step-like enhancement and low-field dip-like 
feature, are the same as those for x = 0.50 crystals. The 
essential similarity of the k(H) behaviors for different 
oxygen contents seems reasonable if the phonon scatter- 
ing by Gd moments is playing the key role, since the oxy- 
gen content apparently has no major impact on the mag- 
netic state of the Gd ions. This result clearly indicates 
the separable effects of free spins and impurities/defects 
on the low-T heat transport of GBCO. It is worth not- 
ing that the quantitative x dependence of k(H) / k(0) is 
rather complicated, namely, the magnitude of k(H)/k(0) 
is slightly enhanced upon increasing x from 0.00 to 0.30, 
while it becomes much more pronounced upon further 
increasing x to 0.50. Obviously, for quantitative descrip- 
tion of the all experimental observations, not only the 
spin-phonon scattering but the oxygen-induced lattice 
disorder must be taken into consideration. 

Although the detailed mechanism of the large mag- 
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netothermal effect in GBCO calls for further investiga- 
tion, one can already capture some important informa- 
tion from the above experimental results. First, the 
phonon heat transport can be strongly dependent on 
the magnetic field even at very low temperature for 
the materials containing paramagnetic moments, like 
high-T c cuprates 12 and other strongly correlated com- 
pounds. Second, if the paramagnetic moments have a 
small enough splitting between the spin states, the spin- 
phonon scattering can survive down to very low temper- 
ature and prevent phonons from entering the boundary 
scattering limit. This makes the data analysis of the low- 
T thermal conductivity very difficult, especially when the 
boundary scattering limit is required for separating the 
electron and phonon terms . 19 i 20 i 21 i 22 



single crystals down to very low temperatures. The main 
features of the magnetothermal conductivity, i.e., a high- 
field enhancement and a low-field dip, can be well under- 
stood in terms of the phonon scattering by the nearly 
free Gd spins. The present finding demonstrates the po- 
tentially significant role of the spin-phonon coupling in 
correlated oxides. 
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